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Linear optics quantum computing (LOQC) is a leading candidate for the implementation of large 
scale quantum computers. Here quantum information is encoded into the quantum states of light 
and computation proceeds via a linear optics network. It is well known that in such schemes there 
are stringent requirements on the spatio-temporal structure of photons - they must be completely 
indistinguishable and of very high purity. We show that in the Boson-sampling model for LOQC 
these conditions may be significantly relaxed. We present evidence that by increasing the size of the 
system we can implement a computationally hard algorithm even if our photons have arbitrarily low 
fidelity and purity. These relaxed conditions make Boson-sampling LOQC within reach of present- 
day technology. 



I. INTRODUCTION 

Linear optics quantum computation (LOQC) [TJ [2] has 
emerged as one of the leading candidates for the imple- 
mentation of scalable quantum computation (QC) [3]. 
Here information is encoded into single photon states, 
and the computation proceeds via a network of linear op- 
tics elements, complemented by measurement and feed- 
forward. It is known that LOQC is universal for quan- 
tum computation [T]. However, owing to fast feedforward, 
the required technology is challenging and well beyond 
the capabilities of present-day experiments. An alternate 
approach, known as 'Boson-sampling', was recently pre- 
sented by Aaronson & Arkhipov (AA) [I], which does 
away with fast feedforward, requiring only single pho- 
ton states, a passive linear network and photo-detection. 
While not believed to be universal for quantum computa- 
tion, AA presented strong evidence that such schemes im- 
plement an algorithm classically hard to simulate, mak- 
ing it of interest to the quantum computing community 
who wish to demonstrate devices with capabilities be- 
yond classical computers. 

Like any QC architecture, LOQC is plagued by difficul- 
ties. In the Boson-sampling model the dominant sources 
of errors are loss, detector and source inefficiency, im- 
purity of photons, and photon distinguishability (caused 
by mode- mismatch or non-identical photon sources) . Re- 
cently Rohde & Ralph [5J considered the issue of loss 
in the Boson-sampling model, and presented evidence 
that within realistically achievable bounds, lossy Boson- 
sampling computers remain classically hard to simu- 
late. In this paper we consider the spectral structure 
of photons, another limiting factor in the implementa- 
tion of Boson-sampling. Conventional wisdom is that 
LOQC places stringent demands on the fidelity and pu- 
rity requirements of photons, which is very technologi- 
cally challenging. We present evidence that in the set- 
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ting of Boson-sampling these tough requirements may 
be relaxed, allowing us to implement a computationally 
hard algorithm even with photons of low fidelity and pu- 
rity. This makes elementary demonstrations of Boson- 
sampling much more foreseeable than other LOQC pro- 
tocols. 

Closely related to Boson-sampling is the quantum walk 
model [5HS], which has recently attracted much experi- 
mental interest [10lll4j . It was argued in Ref. [15 that a 
multi-walker photonic quantum walk is in fact equivalent 
to Boson-sampling. Thus, our results can be interpreted 
as applying to multi-walker quantum walks also. 

The results we present apply only to the restricted 
Boson-sampling model for LOQC and are not applica- 
ble to universal LOQC schemes such as that by Knill, 
Laflamme & Milburn [T], or recent measurement-based 
protocols [T6 1 IT7 ] . 

II. BOSON-SAMPLING 

In the Boson-sampling model we begin with an n pho- 
ton state across m modes, of the form 

IV'in) = |ll,...,l„,0„+i,...,0 TO ), (1) 

where m = 0(n 2 ). The input state passes through a lin- 
ear network, which applies a unitary map to the photon 
creation operators, 

m 

a\^J2 U ^ (2) 

where a] is the photon creation operator associated with 
spatial mode i. In general, in an occupation number rep- 
resentation, the output state is of the form 

\i> out )=J2xs\n[ S) A S K..,nW), (3) 
s 

where YlT=i n< f' > = n ^S, S denotes a photon number 
configuration at the output, and the number of config- 
urations I SI = ( rn+n ^ 1 \ grows exponentially against the 
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number of photons. The probability of measuring some 
configuration is P(S) — \xs\ 2 ^ Per(As) for some ma- 
trix ^5 = f(U,S). Calculating matrix permanents is 
known to reside in the complexity class #P-complete, 
a class strongly believed to be classically hard to solve, 
giving rise to the believed hardness of Boson-sampling. 
See Ref. [1] for a much more rigorous complexity argu- 
ment. The model is illustrated in Fig. [T] We will denote 
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FIG. 1: Model for Boson-sampling linear optics quantum com- 
putation. We begin with some number of modes m. The first 
n modes are initialised with a single photon state and the 
remainder with the vacuum state. The input state passes 
through a passive linear network comprising beamsplitters 
and phase-shifters. We repeat the experiment many times, 
sampling the output probability statistics. 

an instance of a Boson-sampling computer with n pho- 
tons as BosonSamp(n). 

It is known that a universal quantum computer can 
simulate Boson-sampling, but the converse is not be- 
lieved to be the case. Specifically, it is known that 
BosonSampP C SampBQP, but strongly believed 
that BosonSampP C SampBQP. 



III. THE MODE STRUCTURE OF PHOTONS 

In most usual treatments, one represents Fock states 
as \n), or using some power of photon creation opera- 
tors aF In the present study we are interested in the 
temporal/spectral properties of photons, thus this repre- 
sentation no longer suffices. Instead we follow Ref. [18] 
and employ the mode operator formalism for represent- 
ing photonic states. Here we replace the photon creation 
operator with a mode operator defined as 



A ] - 



ijj((j)a]j (w) du, 



(4) 



operators that generate photons with a well defined spec- 
tral structure. Such a spectral decomposition of photonic 
states has been employed previously in a variety of situa- 
tions [TM2"5] . A similar decomposition could be employed 
in the time-domain, simply by taking the Fourier trans- 
form of tp(tj), in the x — y spatial degrees of freedom using 
a double integral, or in the polarisation degree of freedom 
by considering a two-element discrete sum as opposed to 
an integral. But we will focus just on the spectral struc- 
ture of photons for ease of exposition. 

Next, following Ref. [13], we can decompose the spec- 
tral distribution function into a discrete orthonormal ba- 
sis of functions 

^H^^A^M, (6) 

i 

where the coefficients Xi can be calculated from 

A, = J &H>( w )dw. (7) 
For orthonormality, satisfy 

&(w)*&(w) du = (0|%4.|0) = S itj . (8) 
Then, our mode operators can be re-expressed as 

i J 

= ,• W 

i 

Any basis satisfying the constraint from Eq. [8] could 
be employed. Examples could include Fourier bases, 
wavelet bases, Hermite polynomial bases, or bases of fre- 
quency or temporal delta functions. For our study we will 
remain general and not restrict ourselves to any particu- 
lar basis. 

Using such a decomposition one can easily define met- 
rics such as the overlap between two photonic states, 

o M = (oi^tjo) 

tpx (u})*ip2(u) duj 



where A^ j is a creation operator creating a photon with and the fidelity can be defined as 
spectral distribution function tj){ui) in spatial mode j, 
Oj(u>) is a single frequency photonic creation operator 
in spatial mode j, and ^{ui) is a normalised spectral dis- 
tribution function satisfying 



= 10 



(11) 



For spectrally mixed states, in a density operator for- 
malism we write 



|^(w)| 2 dw = 1. 



(5) 



All our integrals are implicitly in the range (0, oo). Thus, 
the mode operators can be regarded as photonic creation 



i 



(12) 
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That is, the state is a mixture over different distributions 
Note that each sub-distribution is pure. The purity 
of a mixed single photon state is 



tr (4,i) =E 



(13) 



IV. 



BOSON-SAMPLING WITH ARBITRARY 
PHOTONS 



We now consider the effects the spectral properties 
of the input photons have on the operation of Boson- 
sampling, considering the cases of both spectrally pure 
and spectrally mixed photons. 



A. Spectrally mixed photons 

We begin by considering the effects of spectrally im- 
pure/mixed photons on the operation of Boson-sampling. 
Our input state can be expressed 



Pu 



)p^j,i 



3=1 



j=l i 



(14) 



Expanding this expression we obtain a mixture of dif- 
ferent combinations of n photon input states distributed 
across the different spectral basis functions. As a simple 
illustrative example, suppose our Boson-sampling com- 
puter has n = 2 photons, m = 2 modes, and the spectral 
distribution functions are supported by two basis func- 
tions £i and £2- Then the expansion of the input state 



from Eq. 14 is of the form 



Pin = 7l,l7l,2 ' P£i,l 

+ 71,172,2 • /0£t,i ® p&a 

+ 72,l7l,2 • P&A ® P£i,2 

+ 72,172,2 ' P&,1 ® P? 2 ,2- 



(15) 



Next, note that because {^} form an orthonormal ba- 
sis, if two photons are present, each in the same spec- 
tral basis state, they will interfere as expected in the 
ideal case, whereas if they are in different spectral basis 
states, they will not. Thus, p^.i <g> p^ 1 ^ gives rise to a 
BosonSamp(2) computer, while p^ lt i <8 P{ 2 ,2 gives rise to 
two independent instances of BosonSamp(l) computers. 
Therefore, Eq. |15| gives us a probabilistic mixture of two 
BosonSamp(2) computers and four BosonSamp(l) com- 
puters. Upon measuring the output statistics, we are clas- 
sically sampling across these multiple instances of quan- 
tum Boson-sampling. 

In general, we may expand Eq. [Ml as 



E IT* 

v£V \j=l 



v ji3 



(16) 



where V is the set of all vectors of length n, with integer 
indices spanning the support of the spectral basis, v is 
one such vector, and v i is the rth element of that vector. 
v can be interpreted as an instance of Boson-sampling, 
and Vi as which spectral basis function the ith photon is 
in for that instance. An illustrative example is shown in 
Fig. [2] 
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FIG. 2: An example instance vector v = (1,2, 0,1) for an 
n — m — 4 Boson-sampling computer. There are three spec- 
tral basis states, the wave-packets shown in grey, each de- 
fined by some mode-function f j, in this case localised temporal 
packets. For a particular v, the photons will occupy some com- 
bination of the allowed spectral basis functions, shown in dark 
grey. The total input state pi n is given by a summation over all 
allowed v £ V, weighted by products of the coefficients jij. 
In this example, v gives rise to two BosonSamp(l) comput- 
ers (in modes £0 and £2 respectively) and one BosonSamp(2) 
computer (in mode £1). Thus, for rihard = 2, we postulate 
that this instance v is hard, which occurs with probability 

p(v) = 71, 172,270, 371.4- 

We are interested in the situation where our Boson- 
sampling computer is classically hard to simulate. To for- 
malise this, let us pick some threshold number of indis- 
tinguishable photons, rihard) upon which Boson-sampling 
is defined to be hard, and we assume rihard < n. Thus 
BosonSamp(n') is hard for n' > rihard- In our expression 
for the spectrally mixed input state we have a mixture 
of many different terms across a basis of different spec- 
tral distributions. Upon measurement we are classically 
sampling across different Boson-sampling problems. For 
classical hardness we desire that the probability we are 
sampling from a classically hard problem be above some 
threshold, p h ard > e. 

Let #(r>) denote the largest number of repetitions 
of a single integer in the vector v. For example, 
#(1, 2, 1, 2) = 2 and #(1, 2, 3, 2, 2) = 3. Then v denotes 
an instance of hard Boson-sampling when #(r>) > rihard- 
In Eq. [16] we are classically sampling across many Boson- 
sampling problems, one for each term in the summa- 
tion. Specifically, the term associated with vector v gives 
rise to a BosonSamp(#:(r))) computer with probability 
p( v ) = Y\.j=ilvj,j, an d some additional smaller Boson- 
sampling computers if #(i>) < n. 

The probability that we sample from a hard instance 



of Boson-sampling is thus given by 

n 

Phard = II 

« I #('")>™hardi = 1 
v I #(f)>«hard 



(17) 



Let us consider the special case where all of the photons 
are spectral mixtures, but identical. In this case jij is in- 
dependent of j. First consider the limiting case where all 
the photons are identical and spectrally pure. In this case 
we can set £ = ip and there will only be one non-zero co- 
efficient in the spectral decomposition, i.e. 7, = Sifi. Then 
p(v) = 1 when v = (0, . . . , 0), otherwise p(v) = 0. If we 
consider the most restrictive case where nhard = n, then 
there is only one surviving term in the summation when 
v = (0, . . . , 0), for which p(v) — 1, and thus Phard = L 
Therefore, in the limiting case of spectrally pure, iden- 
tical photons, our Boson-sampling computer is always 
implementing a hard algorithm, as is expected. 

Next we consider the other limiting case where the 
photons are identical, but maximally spectrally mixed 
across 6 spectral basis states. In this case p(v) = 1/6™ Vw 
and the single photon purity is V = 1/6. A closed-form 
expression for phard is challenging, but a lower bound is 
easily obtained, 



Phard 



> 



E 

fc="hard 



v k (i-vy 



(18) 



For maximally spectrally mixed states, the lower 
bound on Phard is shown in Fig. [3j The important fea- 
ture of this plot is that even for highly impure photons 
we can achieve computational hardness if we use enough 
of them. Thus, for a given degree of hardness there is 
a tradeoff between the single photon purity and the re- 
quired number of input photons. Specifically, as n — > 00, 
the required purity V — > 0. Thus, for identical but mixed 
photons of arbitrarily low purity we can always achieve 
computational hardness for sufficiently large systems. 

The intuition behind this result is straightforward. For 
computational hardness we require that there exists some 
spectral mode within which reside at least nhard photons. 
If the spectral basis states are being populated ran- 
domly, then clearly the probability of this occurring must 
asymptote to unity for large n. 

This is an important observation for experimentalists, 
who are often limited by their photon source technology. 
For example, when preparing single photons by heralded 
parametric down-conversion (PDC), spectral correlations 
between the signal and idler photons will manifest them- 
selves as spectral mixing in the heralded photon. For 
this reason much effort is invested into engineering PDC 
sources with separable spectral distributions. Addition- 
ally, uncertainty in the timing of the pump pulse will lead 
to temporal mixing. However, this result suggests that by 
scaling up the size of the system, such limitations may be 



overcome and computational hardness achieved nonethe- 
less. 



«hard=10 




FIG. 3: (Colour online) Lower bound on phard for maximally 
spectrally mixed, but identical photons. For sufficiently large 
numbers of photons we can achieve computational hardness 
for arbitrarily low purities. 



B. Spectrally pure photons 

Next we turn our attention to spectrally pure photons, 
which in general have distinct spectral structures. The 
input state with arbitrary spectrally pure photons can 
be expressed 



i=i 



j=i i 



(19) 



In the spirit of Eq. [16] we may re-express our input 
state as 

(n n \ 

II A, , A 0), (20) 

3=1 j'=i 3 J 

which is a sum over mutually orthogonal basis states. 
Again we desire that the total probability contributed by 
hard terms is above some threshold, so we define 



Phard 



E 
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E 
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1 I ^Vj,3 
3 = 1 

p(v). 



(21) 



Now we can consider two limiting cases. First, in the 
case of identical photons, applying the same reason- 
ing as before, there is only one surviving term when 
v = (0, . . . , 0), and then phard = 1 as expected. The other 
limiting case is when all the photons are completely dis- 
tinguishable and therefore evolve independently and do 
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not interfere with one another. That is, each resides in a 
different basis state. In this case there is only one term 
in the superposition, corresponding to v = (1, . . . , n), i.e. 
the ith photon resides entirely in £j. Obviously the con- 
dition > Jihard is not satisfied (except in the trivial 
case where nhard — !)• Thus with completely distinguish- 
able photons we are always implementing a classically 
easy computation (n instances of BosonSamp(l)), which 
is also expected. 

We have considered the best- and worst-case scenarios 
for pure photons. We now consider the intermediate case, 
where the photons are pure but have some arbitrary de- 
gree of distinguishability. Deriving a completely general 
expression for this is prohibitive owing to the complicated 
combinatorics. Instead we will derive bounds on the op- 
eration of Boson-sampling where the worst-case photon 
distinguishability is known. That is, the minimum fidelity 
between any pair of photons is known, but we don't know 
the exact spectral decomposition for every photon. This 
is insightful as determining the fidelity between a pair 
of photons is relatively straightforward using simple in- 
terferometors, but performing full tomography of single 
photon states in the spectral degree of freedom is rather 
complex |25j and to our knowledge has never been exper- 
imentally demonstrated. 

Let us introduce a distinguishability parameter 
< a < 1, which captures the overlap between two pho- 
tons' spectral decompositions, a is related to the fidelity 
as J 7 = a 2 . Let .Fmin represent the worst-case fidelity 
between any two photons in our Boson-sampling system. 

For a given value of F r „in, the best-case operation of 
the system is when all the photons are identical, except 
for one which has overlap .Fmin with the remainder. In 
this case the input state is 



\lpin) = (' 



aA 



t 



a 2 A 1 



£1,1 



U A LM- ( 22 ) 



i=1 



Now there are two allowed state vectors, v = (0, 0, . . . , 0) 
and v = (1, 0, . . . , 0), which occur with probabilities 
p(0, 0, . . . , 0) = F mill and p(l, 0, . . . , 0) = 1 - F min . Since 
this represents the best-case scenario, substituting into 



Eq. 21 we have 



Phard < 1, 



(23) 



for n hard < n. 

Alternately, for given F m in) the worst-case operation 
of the system is when all photons have overlap F m i n with 
all other photons. In this case we may write our input 
state as 



n 

U{ aA L + ^*~ 2A L)\°y ( 24 ) 



Since this represents the worst-case scenario, substituting 
into Eq. 121] we obtain 



Phard 



> 



n—k 



(25) 



This lower bound on the hardness probability is illus- 
trated in Fig. |1J demonstrating that high p„ard may be 
achieved even with low worst-case fidelities, provided the 
system is sufficiently large. Again, this is a useful observa- 
tion for experimentalists, who have inherent limitations 
in their photon fidelities. Indeed, as n — >■ oo, the required 

F min ^ 0. 
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FIG. 4: (Colour online) Lower bound on phard against worst- 
case single photon fidelity and number of photons, for spec- 
trally pure photons. For sufficiently large numbers of photons 
we can achieve computational hardness for arbitrarily low fi- 
delities. 



Our combined bound on p n ard is now 



E ( J^min(l-^mi„)"- fc <Pha r d<l, (26) 

for n na rd < n. Fig. [5] illustrates examples of the regions 
satisfying this inequality 
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FIG. 5: Regions satisfying the inequality from Eq. |26| J-min is 
the minimum fidelity between any pair of photons in the sys- 
tem, (left) As our requirement for hardness increases (larger 
Whard), so does the upper bound on the required worst-case 
fidelity, (right) As the number of photons increases, the upper 
bound on the worst-case fidelity required to achieve a given 
degree of hardness decreases. 
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C. Interpretation of hardness 



Importantly, the regions satisfying phard > e in Eq. 17 
and Fig. [3] are not provably hard. We cannot rule out 
the possibility there may exist some 'tricks' enabling ef- 
ficient computation for the mixed states in question (see 
Ref. |5_ for further discussion on this issue). For exam- 
ple, there are known quantum states with large average 
photon number (i.e. n 3> nhard), but whose sampling 
problems are nonetheless classically easy (e.g. coherent 
states and other Gaussian states [26], [27]) . Thus, the re- 
gions within Phard > c m &y be classically hard, while 
the regions outside Phard > e are definitely not classically 
hard. That is, Phard > e is definitely a necessary condition 
for computational hardness, but may not be a sufficient 
condition. 

Similarly, in the case of pure states, Eq. [21] and Fig. [4] 
do not specify provable hardness. Rather, we argue that 
the net contribution from terms which are individually 
hard is above some desired threshold. Of course there is 
in general interference between these different terms, so 
we cannot rule out the possibility that such interference 
effects make classical simulation easier. 

Another important point is that the computational 
hardness of a Boson-sampling device is not just a func- 
tion of the number of input photons n, but also of the 
unitary map U. Even with an ideal input state as per Eq. 
[T] there are some unitary maps U which are always com- 
putationally easy to simulate. The obvious examples are 
permutation matrices, which simply remap input modes 
to output modes without inducing any kind of interfer- 
ence effects. Obviously such a system can be trivially 
simulated for any input. 



D. Overcoming spectral effects 

In universal QC schemes, error correction and fault tol- 
erant protocols can be employed to overcome the effects 
of errors and allow computation to proceed [3]. However, 
in the limited Boson-sampling architecture there are no 
known error correction techniques. Thus fault tolerance 
may not be possible. One technique that is widely used 
in LOQC experiments to overcome mode-mismatch is to 
employ narrowband spectral filtering. That is, we em- 
ploy frequency-resolving detectors (or detectors with a 
frequency filter) and post-select on events where the pho- 
tons are within a narrow frequency range. This has the 
effect of projecting the photons onto a spectral struc- 
ture whereby the different photons appear indistinguish- 
able. While this technique is very useful and widely em- 
ployed in elementary demonstrations, it is very limited 
since post-selection/filtering is equivalent to loss - it dis- 
cards a large part of the wave-packet, thereby reducing 
the detection probability. Thus, if there are n photons in 
the system and the single photon efficiency is p, the prob- 



ability of detecting all n photons is p n , which drops expo- 
nentially against n. In Boson-sampling we ideally wish to 
use a large number of photons to achieve complexity be- 
yond classical capabilities. Thus, due to its unfavourable 
scaling properties, the spectral filtering technique is not a 
satisfactory approach. However, our results suggest that 
if computational hardness is the only objective, then er- 
ror correction may not be necessary. Rather, we simply 
scale our system to have more modes and more photons. 



V. CONCLUSION 

We have considered the operation of the Boson- 
sampling model for linear optics quantum computation 
where the photons have arbitrary spectral structures, 
considering both the cases of spectrally pure and spec- 
trally mixed photons. We derived analytic conditions for 
the relationship between the potential hardness of the 
computation and the spectral structure of the input pho- 
tons. 

We observed that with spectrally impure photons of 
arbitrarily low purity, computational hardness may be 
achieved by scaling up the size of the system. And for 
spectrally pure photons, the worst-case pairwise fidelity 
can be used to construct a lower bound on the hardness 
probability of the computer, and computational hardness 
may also be achieved for sufficiently large systems even 
with arbitrarily low fidelities. 

Our results suggest the otherwise stringent require- 
ments on photon indistinguishability in optical quantum 
computing schemes may be significantly relaxed at the 
expense of a larger system. 

While we have not presented formal hardness proofs, 
our results build on those of Aaronson & Arkhipov, and 
provide circumstantial evidence that Boson-sampling re- 
mains hard even with highly imperfect photon sources, 
provided we scale our systems sufficiently. Therefore, 
present-day limitations in photon engineering technology 
needn't prevent us from constructing devices with capa- 
bilities beyond those of classical computers, and demon- 
strating such devices may be realistic in the medium- 
term. 
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